ULB-TH/12-11 
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In dark matter (DM) models, the production of a y line (or of a "box-shaped" y-ray spectrum) from DM 
annihilation proceeds in general from a loop diagram involving a heavy charged particle. If the charged particle 
in the loop carries also a color charge, this leads inevitably to DM annihilation to gluons, with a naturally larger 
rate. We consider a scenario in which DM candidates annihilate dominantly into gluon pairs, and determine (as 
far as possible, model-independent) constraints from a variety of observables: a) the dark matter relic density, 
b) the production of anti-protons, c) DM direct detection and d) gluon-gluon fusion processes at LHC. Among 
other things, we show that this scenario together with the recent claim for a possible y line from the Galactic 
center in the Fermi-LAT data, leads to a relic abundance of DM that may be naturally close to the cosmological 
observations. 



I. INTRODUCTION 

The trouble with indirect searches for dark matter is that in 
general the possible signatures, like gamma rays and antimat- 
ter in cosmic rays, may not be easily discriminated from those 
of more mundane, astrophysical sources. Two notable excep- 
tions are the possibilities that dark matter may produce signif- 
icant gamma-ray monochromatic lines JTJ, or that it may be 
captured by the Sun or the Earth [3], where it may produce 
high-energy neutrinos. These "smoking guns" for dark mat- 
ter, which have no expected astrophysical counterparts, are 
actively sought for by various experiments but with, so far, no 
conclusive result. In particular the Fermi-LAT collaboration 
has put limits on the possible monochromatic flux [4| from 
dark matter annihilation or decay (see also 0). Recently, it 
has been claimed that the current Fermi-LAT data may actu- 
ally contain some indication of a monochromatic gamma-ray 
signal in the vicinity of the Galactic centre, with Ey ~ 130 
GeV J6). This exciting, albeit tentative result, has been al- 
ready challenged in Q, where it is claimed that the excess 
may have an astrophysical explanation, possibly related to the 
so-called Fermi bubbles regions (see also [8|). Some support 
for a monochromatic line is given in ||9] [10), in which it is 
furthermore claimed that the signal is actually consistent with 
dominant emission from the Galactic centre. This exciting 
possibility has already attracted some attention in the recent 
literature ifTWTl . 

To produce a y-ray spectrum from DM annihilation with sharp 
features, such as a y line (or more generally a narrow box- 
shaped spectrum |[T2l ). in general one needs to invoke a 
Feynman diagram with a charged intermediate particle (see 
e.g. [22-24| for non-SUSY instances), which can manifest it- 
self only at the one-loop level (an exception to this rule of 
thumb is the scenario of Ref. ll25ll ). If we limit ourselves to 
one-loop induced couplings, then there are two topologies for 
a y line, see FigJTJ either through a quartic vertex, or through 
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an s-channel annihilation into an intermediate particle S that 
couples to two photons. The former is an effective operator 
of dimension 7 (6) for a fermionic (scalar) DM candidate. 
While the two processes become indistinguishable in the limit 
ms 3> 2moM, the quartic coupling does not necessarily require 
a UV completion with an intermediate particle: a box diagram 
for instance, like in fl4l . or in the MSSM Il26ll27l . is another 
common possibility. 




FIG. 1 : Topologies with pairs of gamma-ray lines. The black dot is 
meant to represent a one-loop induced coupling. 

Similarly for a box-shaped spectrum, were the particle(s) that 
produce the y are on-shell, there are three topologies (Fig|2]) 
still at one-loop: t, s or quartic channel annihilation of DM to 
two on-shell particles, or messengers, that we call generically 
S (there may be more than one S), with subsequent one-loop 
two-body decay of S. The box-shaped scenario, unlike the y 
line scenario, has the particularity that the production of hard 
y is not necessarily one-loop suppressed, as it is proportional 
to the tree-level production of on-shell S particles, and to their 
branching ratio into photons (this is a particular instance of 
secluded DM B4l ). This allows to produce a large flux of 
y, even if the couplings are small. Familiar examples of S 
neutral particle decaying to photons are scalar fields (tc , ax- 
ions, Higgs) which decay into y through a loop diagram with 
charged particles. 

In the following we analyse the consequences of a scenario 
where some of the charged particles in the loop carry color, 
in which case one has inevitably similar diagrams with a pair 
of gluons instead of photons. One important property of these 
one-loop diagrams is that the ratio of gluon-to-photon pro- 
duction has little dependence on the topology one considers 
in Fig. 1. Another one is that the production of gluons is nat- 
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FIG. 2: Topologies giving to a box-shaped gamma spectrum. Anni- 
hilation proceeds through the decay of a, on-mass-shell, messenger 
particle. A larger coupling to gluons implies that the dominant con- 
tribution comes from decay of the messengers into two gammas pairs 
and two gluons pairs. 



urally larger than the production of photons. As an illustra- 
tion, let us consider the example of coupling to gluons and 
photons pairs through a one-loop, triangular diagram with a 
single heavy vector quark of electric charge Q e , and no weak 
isospin. In this case the ratio of gluons to prompt gammas 
produced per annihilation is 
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with C a numerical constant that, in the heavy quark limit 
(m Q > m DM ), is 1 
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For the numerical estimate we have set Q e = 1 and have we 
taken as and a at the scale moM = 130 GeV. 
In the second equality of ([T]) we have not taken into account 
the contribution from the (ov)yz and (ov)zz channels. These 
are automatically present, but they are sub-leading and we ne- 
glect them for the sake of our argument. Now, if there are no 
other annihilation channels than the one-loop processes dis- 
cussed so far, then the total annihilation cross section may be 
approximated by 
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so that, unless the electric charge is larger than Q e = 1, one 
gets naturally a large ratio. 

Similarly, in the case of a narrow box-shaped spectrum, one 
gets (neglecting again the sub-leading decay of S to yZ and 
ZZ): 
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just as in ([T]). Taking into account that C is naturally much 
larger than one, most of the photons are produced in the pro- 
cess DM DM— > gg + yy, which implies that 
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The results of Eq. (|3]l and Eq. |5]l are interesting in two ways. 
First of all, one can compare these ratios with the ratio of the 
canonical cross section required for thermal freeze-out, 
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to the yy cross section for the 130 GeV line claimed in Ref.|6] 
(for an Einasto profile), 
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This ratio equals ~ 25, very close to the values quoted above 
for Q e — 1 ! Therefore, the scenario we consider may give a 
natural explanation for the strength of such a y line. Second, 
and this is of broader scope, it is important to confront such 
a scenario to other indirect, direct and collider constraints. In 
particular, the existence of a large channel into gluons may 
give rise to a plethora of other particles, like anti-protons, and 
also opens interesting possibilities for direct detection of dark 
matter, as well as colliders signatures. 

In the sequel, we first derive model independent constraints 
from the observed anti-protons flux (section|II]). Then we con- 
sider more model dependent constraints that may be set by 



direct detection experiments, section III Finally we consider 
colliders constraints (jets+missing energy) and possible sig- 
natures in section UV] 



II. ANTI-PROTON FLUX 

For DM particle with mass above the GeV scale the produc- 
tion of anti-protons from gluons is unavoidable in the scenario 
we consider. Then, if we take for granted a positive gamma- 
ray line signal, an anti-proton flux is the most straightforward 
outcome, as the only relevant parameter is the ratio of Eq.([T]i. 
More generally, given a scenario with both gamma and gluon 
lines, with a specific prediction for the ratio C, we may ask 
the following question: given the maximum anti-proton flux 
allowed, what is the prospect for gamma-ray line detection? 
In both cases, the constraints are independent of the topology 
considered. Also, independently of a y line, we may ask what 
is the maximum annihilation cross section into gluons allowed 
by the anti-proton data, and for which DM mass is it equal to 

(CTV) F0 ? 

In this section we refer specifically to bounds on annihilation 
into gluon and gamma-ray lines, hence to Eq.|3]l 2 . In Figj3]we 



Notice that, if the 0(a.s) corrections to yy production are small, there are 
substantial for decay into gluons, multiplying the ratio of Eq. J5J by a factor 
of about 2 (see e.g. |43 ]). Unless stated, we do not consider explicitly these 
NLO corrections, as our approach is essentially model independent. In 
other words, such corrections are included in the definition of C. 



2 For secluded dark matter, each annihilation produces four photons with 
energy equal to moM/2. Hence, in Figs.Q and |5J, the horizontal line 
that corresponds to freeze-out abundance should be replaced by (o~v) = 
{o"v)fo/2, and the values on the horizontal axis should be multiplied by 2. 
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give, for illustration, the flux of anti -proton obtained assuming 
a single heavy vector quark of charge Q e = 1 and saturating 
the flux of gamma to that of the gamma line of Ref.[6 1 at 130 
GeV. In this plot, we have taken into account the contributions 
from both the yy and yZ vertices. For the calculations of the 
flux of anti-protons we have used a popular semi-analytical 
approach, as summarized in Ref. |28l , to which we refer for 
more details. The figure displays two possible fluxes, using 
the standard MIN, MED and MAX sets of cosmic ray propa- 
gation parameters and NFW and Einasto profiles for the distri- 
bution of dark matter in the Galaxy. The data points are those 
collected by the PAMELA collaboration (29). Basically, the 
conclusion is that the MIN and MED curves are consistent 
with the data, while the flux corresponding to the MAX pa- 
rameters overshoots the data and so is excluded. 3 
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FIG. 3: Anti-proton fluxes from a WIMP candidate of iudm =130 
GeV and annihilating in gluon pairs, compared to the PAMELA mea- 
surements. The fluxes are given for the NFW and Einasto profiles 
and for MIN, MED and MAX cosmic propagation parameters. 

Independently of the 130 GeV gamma-ray line, we give in 
Fig |4] the bounds on annihilation into gluon lines allowed by 
the PAMELA anti-proton flux, taking into account the sec- 
ondary anti-protons background flux produced by cosmic ray 
interactions |30|. Depending on the astrophysical setup, the 
bounds on annihilation cross-sections in gluon lines may be 
larger (or smaller) than the canonical freeze-out value. The 
standard way to read this plot is that, depending on the prop- 
agation parameters, a WIMP candidate annihilating domi- 



3 A priori, another possible signal are positrons in cosmic rays, but we have 
checked that the flux of positrons is always substantially smaller than the 
PAMELA data. We have also checked that the overall flux of photons, 
including those from gluons, is consistent with various indirect constraints 
(diffuse galactic and extra-galactif gammas fluxes). 



nantly into gluon lines, i.e. (cv) gg ~ (gv)fo, is excluded be- 
low a given mass: moM ~ 180 GeV for MAX, moM ~ 80 
GeV for MED, m DM ~ 16 GeV for MIN (not shown). One 
should however keep in mind that a larger cross section may 
be implemented in the Galaxy (at the prize of some fine-tuning 
and/or more model building) due to an astrophysical boost 
factor (however numerical simulations limits these to be of 
order a few [31 1), or a particle physics effect (Sommerfeld en- 
hancement, resonance,...). Our attitude here is agnostic, and 
we take these bounds as the maximum possible values allowed 
by current data (and the modelling of the CR background). 
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FIG. 4: Upper bounds on the annihilation cross section of DM into 
gluon pairs as set by the PAMELA anti-protons flux measurements 
for Einasto and NFW profiles, and three set of propagation parame- 
ters. The horizontal line corresponds to the standard thermal freeze- 
out annihilation cross section. 

Now these bounds may be mapped on bounds on the max- 
imum possible gamma-ray line fluxes. Although the map- 
ping is straightforward, we find it useful to illustrate it in a 
separate plot. In Fig|5]we show, for four possible values of 
C = 5,25, 100 and 500, the corresponding maximum annihi- 
lation cross sections in gamma ray lines from Fermi-LAT (at 
95% C.L.). In the same plot we show the estimation from 
Ref. ll32l of the potential reach of the HESS experiment. Not 
surprisingly, a very low value, C = 5, is basically excluded 
for a mass up to the value which has been probed by the 
Fermi-LAT experiment, moM — 200 GeV, even for the MAX 
parameters. At the other extreme with large values of C (as 
exemplified for C = 500), the expected flux are too small to 
be probed by Fermi-LAT. The most interesting situation is the 
one of intermediate values, C = 25 — 100. For instance, for 
C = 100 and MED propagation parameters, a WIMP, which 
emits gluons on top of photons with such a factor C, is al- 
lowed by gamma lines for all masses. The same situation 
holds for MAX parameters and C = 25. Alternatively adopt- 
ing the MED parameters and C ~ 25, Fermi-LAT should have 
seen already a signal with moderate significance, and the 130 
GeV tentative line corresponds to such a case along our sce- 
nario. Ultimately if a y line turns out to be accurately observed 
in the future, the results of Fig. [5] will allow to set an upper 
limit on C (imposing that the annihilation cross section to glu- 
ons does not exceed the thermal freeze-out one). This will 
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FIG. 5: From the results of Fig. 4, upper bounds on the DM annihilation cross section into monochromatic photons, for C = 5, 25, 100, 500. 
The solid (dashed) black line gives the experimental upper bounds from Fermi-LAT assuming a NFW (Einasto) profile |4|. The green line is 
an estimate of the possible reach of the HESS experiment (taken from 1321 ). The horizontal line corresponds again to the FO annihilation cross 
section, rescaled by C. 



provide information on the "beyond the SM" heavy charged 
particles at the origin of this line. The discussion of this sec- 
tion also applies to narrow "box-shaped" spectrum by simply 
rescaling the mass of DM and the ratio C, according to relation 

0- 



III. DIRECT DETECTION 

A potentially interesting feature of the gluon line scenario is 
that it opens the possibility for direct detection through scat- 
tering on nuclei. The connection between the gluon cross sec- 
tion and possible direct detection is not as model independent 
as the anti-protons signature, but depends on the topology and 
on the DM particle nature. For some topologies the relation 
remains one-to-one, with, as we show, prediction of a signal 
of the order the present experimental sensitivity. 

Case 1. y-line, quartic interactions. The most straightforward 
case is that of an effective quartic coupling between two DM 
particles and two gluons, which gives a one-to-one relation 
between the gluon cross section and the direct detection rate. 
In this case, a single scale, A, controls both the annihilation 
cross section and, for the gluon vertex, direct detection and 
signatures at hadron colliders, which we will discuss in sec- 



tion|IV] see Refs. ||33ti37l . 

For our argument, the most relevant operators involving two 
DM particles and two gluons are 
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(9) 



for a scalar (S) and a fermion (F) dark matter candidate respec- 
tively (we consider complex scalar and Dirac fermion DM — 
the generalization to real fields is straightforward). Generally 
speaking, one could also consider fermionic operators involv- 
ing %fXf, but their annihilation is P-wave suppressed, and so 
they are irrelevant (baring a very large boost factor) for in- 
direct detection, and a forteriori for a gamma ray lines sig- 
nal. Similarly, one could think of operators mixing scalar and 
pseudo-scalar quantities, like 
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These operators are allowed but, as they break CP symme- 
try, they are in our opinion somewhat less likely to have large 
coefficients, and for conciseness, we simply discard them. 
Hence we are left with the two possibilities of Eqs.(8]9 1. 



5 



These operators lead to the following annihilation cross sec- 
tions into gluon pairs: 



(av) gg = a 2 s m 2 Xs /(9K 3 A 4 s ) 
(av) gg = a 2 m 4 , F /(2n 3 A 6 F ) 



(10) 
(11) 

Imposing a thermal freezeout abundance, Eq. (If)}, they give 
respectively 



Xs 



130 



134GeV, 

,v ^ 

171GeV 
were we set as = as(130GeV). 



GeV, 



"If 



130 



GeV, 



(12) 
(13) 



<crv>„ ; 




Xenon I T Reach 



FIG. 6: Bounds from XenonlOO (solid) on a scalar candidate (masses 
in GeV) with coupling to gluons through the operator of Eq.jSJ, 
(annihilation cross section in cm 3 s _1 ), the Tevatron limits (solid) 
from jets + missing energy, and the prospects for XenonlT (dashed, 
adapted from |42|) and the LHC (dashed, adapted from 1361 ). For 
av ~ 1 pbarn (horizontal line), candidates heavier than ~ 130 GeV 
(vertical line) are allowed by current data. 

What are the prospects for direct detection (DD)? The op- 
erator Q is a priori out of reach of direct detection experi- 
ments because its scattering cross section is both recoil mo- 
mentum suppressed, and Spin-Dependent (SD). The other op- 
erator however, Eq.([8]l, leads to a standard Spin Independent 
(SI) elastic cross section on nuclei, and so may be tested by 
DD experiments. The procedure to derive the effective cou- 
pling of %s to a nucleon is pretty straightforward [38 1. The 
matrix element of the gluon operator on a nucleon state may 
be written as 
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12% ' 



(N\G"G a » v \N) =f N m N (N\y N v N \N)(=f N m N ), (14) 



with an effective coupling — 2/27 « 0.074. This is pre- 
cisely the contribution that would be induced by a heavy quark 



loop. In our comparison with experiments, we use a more pre- 
cise expression that takes into account the contribution from 
the light quarks to the nucleon mass (see for instance [39]), 

2 

/jv = 27 C 1 - (l+ry/2)a,w/m N ) 

where y is the fractional strange-quark content of the nucleon, 
r = 2m s j (m u + m^) and (%v is the pion-nucleon sigma term. 
There are substantial uncertainties on fy. Our estimate is fy = 
0.053 ±0.011 (at l-o). Regardless, from the above we get 
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where /j is the %s/N reduced mass. 

In Fig|6j we use Eqs.([T0)> and ( 15 1 to constrain the annihila- 
tion cross section into gluons using the exclusion limits set by 
XenonlOO [40]. The horizontal band (reddish) corresponds to 
the standard freeze-out annihilation cross section ~ 1 pbarn 
(see BP for a recent re-appraisal of (<3v)po)- The other band 
corresponds to the exclusion limits set by XenonlOO. The 
width of the band corresponds to a 2-o variation of the pa- 
rameter fy. The vertical line corresponds to m^M = 130 GeV. 
Also, we have allowed, in the relation between the annihila- 
tion cross section and the elastic SI scattering cross section, 
for the possibility of next-to-leading corrections in as to the 
annihilation cross section, which may amount to a factor < 2, 
based on the expression for Higgs decay in gluon-gluon (see 
for instance [43 1). The (corresponding) corrections to the elas- 
tic cross section are more difficult to assess, as the expression 
is non-perturbative in nature. To be concrete, we conclude that 
the XenonlOO data exclude all the candidates below ~ 130 
GeV (except for light, ~ few GeV, candidates, which we do 
not address here): hence a candidate with moM = 130 GeV 
will be tested by the forthcoming XenonlOO data. The future 
XenonlT experiment may probe candidates up to ~ 385 GeV. 



Case 2. y-line, s channel exchange. The operators of Eqs 
and (|9J can be induced by the s-channel exchange of, respec 
tively, a scalar S or a pseudo-scalar P particle (messenger in 
the sequel), Fig{T] For instance the following Yukawa and 
one-loop effective interactions 

L D VsSxhs+yQ^^G^, (16) 
P a? 

L D iyePytFlslF+yQ—^G^G^ , (17) 

will induce such interactions. In the effective coupling be- 
tween the messenger and gluons, we have in mind the effec- 
tive coupling of a single heavy vector quark. More generally 
speaking, mg/yg can be considered as being an effective mass 
scale, just like in operators (8]9 1, which encompasses the con- 
tributions of many degrees of freedom running in the loop. 
Accordingly, in the sequel we will often set yg = 1 . The dis- 
cussion of the relic abundance and of indirect signals depends 
on the mass of the s-channel messenger. 

Clearly, for a heavy messenger, m 2 SP >> 4m^ M , the dis- 
cussion reduces to the previous case, Fig. [6] One recovers 
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Since moM/niQ and 4mf w /m|p are smaller than one, candi- 
dates with moM of order 130 GeV for example, require values 
of the couplings and of /j$ /mDM that are much larger than one, 
typically ~ 10. This condition is also required for the colour 
state Q to be heavy enough to escape to the present collider 
bounds. Although these constraints could be slightly relaxed 
if there were several colored states, this discussion illustrates 
how demanding it is to get a 130 GeV line from the operators 
of Eqs. (|8]l-(|9} when rephrased in terms of an explicit model. 

For the opposite case of a light messenger, 4m 2 JM > m| P , one 
gets Eqs. ( 18 l-(21 1 provided one substitutes m^ p by 4m|, M . 
In this casemeTast suppression factor of Eqs. pO|-pT} dis- 
appears, but still the condition (ov) gg ~ (ov) po requires quite 
large couplings. The main difference with a heavy messen- 
ger is that the elastic scattering on nuclei is enhanced by a 
factor (4m| )M //n| i ,) 2 . Therefore, given the bounds shown in 
Fig. do*|, such a scenario with a candidate of m% s = 130 GeV 
and thermal relic abundance is excluded. Interestingly, for 
the case of a fermionic candidate, if mp is much lighter than 
m Xf , one could get a large enough enhancement to have an 
observable spin-dependent signal, in spite of the q 4 recoil mo- 
mentum suppression. This enhancement is maximum for mp 
of order of the ~ keV recoil energy probed in direct detection 
experiments. Concretely a ratio of (mp/m% F ) ~ 10~ 4 , which 
corresponds to mp ~ 10 MeV for m% F ~ 130 GeV would be 
required to be within reach of current SD direct detection ex- 
periments. 

Finally if m$p ~ 2mp)M, we have two possibilities: 
Firstly, we may have resonant annihilation in the GC. While 
the bounds from direct detection (for a scalar candidate) are 
unchanged, this resonance boosts the signals in both gammas 
and gluons. These are constrained by the bounds on the flux 
of anti-protons and gamma lines, see Figs.(4j5i. If annihila- 
tion in gluons fixes the relic abundance, then, for instance, 
the boost from the resonance is constrained to be less than 
~ 5 (MED) for mp,M = 130 GeV. A corresponding bound on 
the gamma ray line flux may be read in Pig . (j5j : for the same 
setup, C must be larger than ~ 25. For example an 0(100) 
resonance boost of the y signal at 130 GeV would require a 
value of C of order ~ 2500. This is not inconceivable: for 
instance, the heavy fermion may carry color but have a small 



electric charge (take for example Q e = 1/3 in Eq. (|2|); if there 
are more particles in the loops, then there is more freedom, in- 
cluding the possibility of playing with different mass scales. 
Secondly, the resonance may be responsible for the relic abun- 
dance. Then both the direct and indirect signals are depleted, 
and the constraints from direct detection (again for the scalar 
DM) may be much weaker than those shown in Fig. |6]). In this 
case, if we focus on the 130 GeV line, and the relic abundance 
is dominantly from annihilation into gluons, then C must be 
smaller than 0(25), for otherwise (ov) gg > (ov)fo- A way 
to accommodate a low C value is to invoke, not only heavy 
vector quarks, but also heavy vector charged leptons (possibly 
with smaller masses) 4 . Of course, the limiting possibility is 
a scenario with no heavy quarks (and perhaps the relic abun- 
dance from resonant annihilation in y,Z), in which case there 
are no constraints from direct detection. Such a scenario con- 
stitutes of course the extreme counterexample to the perspec- 
tive we follow here. Still, in general one expects that C> 1 
in beyond the SM scenarios. 

Case 3. Box-shaped spectrum, t-channel. In case of secluded 
dark matter, the dark matter relic abundance and the indirect 
signal are not directly related. Consequently, if the discus- 
sion is the same as for y line for indirect signals, the prospect 
for direct detection crucially depends on the S particle life- 
time. Also, here we only consider annihilation of a scalar par- 
ticle, as annihilation of a fermionic candidate in the t-channel 
is P wave suppressed: this has little impact on the relic abun- 
dance, or direct detection, but it precludes indirect detection 
and so is of little direct interest for our purpose. 
In the case of t-channel, the annihilation cross section of a 
scalar DM is given by 
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where our conventions for the coupling between the scalar 
DM and the messenger S are the same as in Eq. (16 1. For 
reference, we give the decay rate of the S particle in yy, yZ 
and gluons for the case of one heavy vector quark of charge 

Q e = h 
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yy) (24) 
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4 For a fixed value of (o"v) gg < (ov) fo, the constraints on the candidate from 
direct detection may be read (again for scalar DM) from Fig.|6j. 
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The spectrum of photon (or gluons) produced per decay of the 
messenger has the shape of a box fl2l . 
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(27) 



For this process, the WMAP relic abundance is obtained for 5 

^- 16 °(26^v) 3/2GeV (28) 

for moM 3> ms- If, on the contrary, the DM and the messenger 
are nearly degenerate, the box spectrum is narrow (i.e. close 
to a line) and a good approximation for (is, valid within 10%, 
is 



^^260^) 



3/2 /OA 



1/8 



GeV, (29) 



where 8 = (ra DM - m s )/m DM . 

The direct detection cross section is that of Eq.( 15 I, with the 



scale A given by Eq.( 18 I. A convenient way to write this cross 
section is 



a si w 3.7- 10 
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/260GeVy ( (i s 



7 



'Is 



T{S- 



'Xs 



(30) 

which is valid in the approximation ms —> m Xs . The ap- 
pearance of the decay rate of the messenger (here in gluon 
pairs) is just a way to re-express the dependence in yg/mg. 
Since ms « m Xs , and (is may be fixed by the relic abun- 



dance, Eq.(29i, the cross section essentially depends on the 



DM mass, and the decay rate of the messenger (and also on 



the mass splitting, through Eq.(29 1). To fix the ideas, we con- 
sider 3 instances: 

1) Take ra Zs . = 260 GeV and (is = 180 GeV (this corresponds 
to e ss 0.01, ms = 257 GeV). For this DM mass, the current 
exclusion limit set by XenonlOO is Osi < 2.2 ■ 10~ 44 cm 2 . This 
translates into an upper bound on T(S — > gg) < 0.32 GeV. 

2) Take m Xs = 60 GeV and (i s = 20 GeV (e « 0.17, m s = 50 
GeV). This candidate is interesting because it is for such mass 
that the XenonlOO constraints are the strongest. Taking the 



current XenonlOO limit a si < 6.8 • 10 45 cm 2 we get that the 
decay rate has to satisfy F(S — > gg) < 1 .5 • 10~ 5 GeV. 
3) Take m Xs = 130 GeV and (t s = 57 GeV (e « 0.04, m s = 125 
GeV). This candidate is relevant for collider constraints (see 
next section). The XenonlOO limit is Osi < 1 .2 • 10~ 44 cm 2 so 
that r(S gg) < 3.2 • 10~ 3 GeV. 

Clearly, the bound from XenonlOO on the m Xs = 260 GeV 
is not very strong (for Jq = 1 it would imply that mg > 1 8 
GeV) 6 . For the lighter candidate, thanks both to the stronger 
limit set by XenonlOO and the scaling of the cross section 
with m Xs and the decay rate in ms and rag, the bounds on the 
scale rag are quite relevant. For yg = 1, m Xs — 60 GeV gives 
ra e > 300 GeV, while for m Xs = 130 GeV, ra e > 60 GeV. 
For reference, the future Xenon IT experiment will bring these 
limits to: 

1) T(S -> gg) < 1.0 • 10~ 3 GeV for m DM = 260 GeV, corre- 
sponding to rag > 300 GeV. 

2) T(S ->gg)<l- 10~ 8 GeV for m DM = 60 GeV, correspond- 
ing to rag > 4.3 TeV The latter constraint is probably stronger 
than the LHC reach. 

3) T(5 gg) < 1.2- 10- 5 GeV for m DM = 130 GeV, corre- 
sponding to rag > 1 TeV. This is similar to the current LHC 
reach, see next section. 

Case 4. Box-shaped spectrum, quartic coupling and s-channel. 
We discuss these two cases only briefly, as the prospect for 
direct detection is either null, or it reduces to one of the 
previous cases. Indeed, if the messenger has a sizable 
vacuum expectation value, then we are back to a gamma-ray 
line through an s-channel (Case 2). If instead the vacuum 
expectation value is small or nonexistent, direct detection 
is doomed by the necessity of exchanging two messengers 
with the nucleon in an elastic scattering process, a process 
that is loop suppressed. A naive estimate gives (for the scalar 
candidate) 



Osi ' 



,3.10-»f^Vf^V 
V m s J \ rag J 



. cm z (31) 
rag J 



where we have only included the dominant parameter depen- 
dence (and have set muM ~ 260 GeV). This is true for the 
quartic coupling topology as well as for the s-channel one, 
with in the later case, perhaps, the possibility of further sup- 
pression of the direct detection signal if a resonance occurs, as 
above. In principle one may have an inelastic process, with, 
say emission of a pair of photons, but we estimate this to be 
below the reach of direct detection experiments. 



5 For the above processes to be relevant, we have to check that the annihi- 
lation in the s-channel of Fig[T] (already discussed in the previous section) 

is sub-dominant. This will be the case provided n$>Ji ( j GeV. A re- 
lated issue with DM annihilating through a scalar messenger is that mixing 
with the Higgs must be small for the dominant channel to be decay into 
gluon pairs. Concretely, for a mj ~ 260 GeV messenger, we estimate that 
the bound on the S-Higgs mixing angle must be 9 < 1()~ 3 \ mQ j' ^ s 
mixing is not prevented by any symmetry, such a small value is not natu- 
ral. A pseudo-scalar would be better, as in this case one could invoke CP 
to prevent mixing with the Higgs. We leave such considerations for future 
work. 



IV. COLLIDER CONSTRAINTS 

The predictions of our scenario at colliders depend, as for di- 
rect detection, on the specific topology of the relevant process. 



' This quantity should be considered as an effective mass scale, not neces- 
sarily the mass of a heavy quark. 
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Case 1. y-line quartic interactions. Again this is the most 
straightforward case, as the interaction with gluons and/or 
photons depends only on one scale parameter. As in section 
III we consider the two operators of Eqs.( 8]9 1, which are the 
most relevant ones for indirect detection. Collider constraints 
on such operators rest on the possibility of creating DM pairs 
in collisions, the most obvious signature being the production 
of jets + missing energy. Taking into account the Standard 
Model background, limits may be set on the scale parameters, 
for each effective operator ll33T[37l . Concretely, we refer to 
the analysis of Ref.|j3"6fl, where the collider constraints have 
been studied. 7 In Figrol we have reported the bounds from 
ll36l . both from the Tevatron and the prospect for the LHC 
(i.e. for Ecm = 14 TeV, and 30 fbarn -1 ). The bottom line is 
that the current collider limits are less constraining (but are in- 
dependent) from those from direct detection searches. 8 On the 
long run, the LHC may exclude candidates < 100 GeV, 
but these will be superseded by Xenon IT (if no DM events are 
observed). Similar, but slightly stronger conclusions may be 
reached for the fermionic candidate, see Fig. (11) in Ref. l36l 
for which, we recall, there are no limits from direct detec- 
tion: the Tevatron data exclude a fermionic candidate below 
m %F < 40 GeV ; the prospect for the LHC is m Xf < 700 GeV, 
assuming the standard freeze-out abundance. 
Case 2. y-line, s-channel exchange. Regarding the jets + 
missing energy constraints, the discussion for the second di- 
agram of Fig. 1 reduces to the previous case if the mass of 
the s-channel messenger particle is heavier than the other en- 
ergy scales, m 2 sp ^> 4m 2 3M and p ^> s where s refers to 
the momentum square of the particle P at LHC. If instead 
m| p is light enough to be produced at Tevatron or LHC, with 

4m DM < m s P' ^ s production is resonantly enhanced which 
strengthens the collider sensitivity (whereas the annihilation 
cross section is not enhanced). If on the contrary 4m^, M ~ mp, 
the annihilation cross section is resonantly enhanced, in which 
case the collider constraints are weaker than those for the 
quartic interaction (Case 1). 

This brings another possible signature at colliders, which is 
the production at the LHC of an on-shell messenger particle 
through gluon-gluon fusion, followed by its decay into a pair 
of gammas. We discuss the prospect for this signal in Case 3, 
for which such a signal is a natural prediction. 
Case 3. Box-shaped spectrum, t-channel. Regarding the jets 
+ missing energy constraints, the same comments as in Case 
2 apply, so we do not repeat them here. 
If we focus on the possibility of a narrow box-shaped spec- 
trum, which implies that mj)M ~ nis,p, then one prediction of 
this scenario is the on-shell production of the messenger par- 



7 The operators of Eqs. j8|9| correspond respectively to the operators C$ and 
D 14 in (36). 

8 The limits from the recent LHC runs are not available for the effective oper- 
ators of DM into gluons, Eqs. j8|9) . The current LHC constraints on other 
operators, like effective couplings to quark-anti-quark pairs, give bounds 
that are only a factor of a few stronger than those based on Tevatron data, 
and weaker than the current XenonlOO bounds (for Eq.(|8j). The same con- 
clusion, we expect, should hold for the gluonic operators. 



ticle at colliders. The process is precisely similar to the pro- 
duction of the Higgs, so we will heavily rest on what is known 
for this particle. The smoking signal is the production of di- 
photons, with invariant mass ~ moM (a clean signal). The 
production of a pair of gluon jets is also potentially interest- 
ing, but the background is much larger. 
The gluon-gluon to gamma-gamma ratio of production cross 
section is directly related to the ratio of Eq.Q. For the sake of 
illustration, we will consider two DM candidates with moM — 
260 GeV and mdm =130 GeV. As before, we consider the 
heavy fermion limit (4mg > mf), in which case the following 
relation is a good approximation, 



cjppjgg) 

<3{PP{gg) - 



>S+X) 
H + X) 



GwKgg 
<5w(gg 



H) 



TeVV 



: 0.045 \y Q - 

V m Q J 

(32) 

where o(pp(gg) H+X) is the Higgs boson production 
cross section through gluon-gluon fusion. Taking the latter 
for m H = 260 GeV, we get from Ref.P31 

c{pp{gg) -^H+X) = 3.12 (13.43) Pb at y/s=7 (14) TeV, 

(33) 

from which we infer that the production cross section of a 
m$ ~ 260 GeV particle S is 

/ TeV \ 2 

°(pp(gg) -+S+X) = 0.14(0.61) [y Q pb (34) 

V m Q J 

at y^s = 7 (14)TeV. The branching ratios of Higgs in gluon 
pairs or di -photons at m# = 260 GeV are small R31 : 



Br(H->rf) ~ 2x 10~ 5 , 
Br{H^gg) pb 7x 10 



(35) 



-4 



while the S decays dominantly into gluon pairs. (We note in 
passing that the ratio of these branchings is Ch ~ 32.) From 
the above numbers, we infer that the number of di-photons 
events from S, normalized to that from a Higgs at m# = 260 
GeV, is 



N(H -> yy) 



2-10 3 / TeVV 

~ — — x [yQ — 

m s =260GeV C V m Q J 



(36) 



We did not find constraints on di-photons in the literature (al- 
though they must exist) for such a large invariant mass scale. 
For comparison, let us consider the sensitivity one may have 
for m$ ~ 125 GeV (corresponding to moM ~ 130 GeV, see 



Section IIIi. In this case we may refer to the current con- 
straints on di-photons from Higgs searches at the LHC 11461 . 
For niH = 125 GeV, we have 



7V(5->yy) 



N(H -> yy) 



m s =l25GeV 



30 / TeV Y 

x I y Q 

V m Q J 



C 



(37) 



An excess of di-photons at such energy is quite constrained 
by the current LHC data, in particular if we assume that the 
Higgs mass is indeed 125 GeV (see for instance Ref. flTl ). If 
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we assume that the ratio of Eq.(37 1 is constrained to be < 1 
by LHC data, we have rag > 1 lev for jq = 1 and C ~ 30. 
This is stronger than the current Xenon 100 bound (see Sec- 
tion [Hi »: rag > 60 GeV (for the same Yukawa coupling) and 



comparable to our estimate of the XenonlT sensitivity reach, 
rag > 1 TeV. We tentatively conclude that the constraints on 
such processes are in a ballpark at the interface of direct de- 
tection and colliders, and are worth being investigated more 
systematically, something we leave for future work. 

Case 4. Box-shaped spectrum, quartic coupling and s-channel. 
If the messenger has a sizable vacuum expectation value, then 
the situation regarding the constraints from jets + missing 
energy is analogous to that of Cases 2 and 3. 
If not, these topologies imply that jets + missing energy are ac- 
companied by the emission of gluons (or photon) pairs from 
the second messenger particle in the diagram. Assuming that 
the latter is off-shell (we do not consider the possibility of 
two different messenger mass scales) we expect that these pro- 
cesses are much suppressed compared to the signals discussed 
in Cases 2 and 3, and so are not very constraining. 
Regardless, we may also create the messenger particle and 
observe its decay into di-photons, as discussed in the previous 
section (Case 3). 



V. SUMMARY AND PROSPECTS 

The observation of a sharp y feature at the level of sensitivity 
of the Fermi-LAT (or future Cerenkov telescope arrays) could 
be explained from annihilation of dark matter particles. Such 
a process is expected to occur through loop diagrams, and 
may be associated to various topologies, which at one-loop 
are shown in Figs |l|2| Such a signal would obviously point 
to a relatively large annihilation cross section into y, which 
could be explained either through large couplings and/or the 
proximity of a resonance (gamma-line) or, without the need 
of any of such features, through a secluded dark matter can- 
didate (box-shaped spectrum). A natural question one may 
ask is: what if the heavy charged particles in the loop, which 
emit the photons, also carry color charge? After all this is the 
case for the classic instances: jiq, Higgs and axions. Annihi- 
lation in gluon pairs is allowed in many dark matter models 
(in particular the MSSM), but suppressed. Here we consider 
the simple possibility that the annihilation is mostly into glu- 
ons. This may arise in a framework in which the dark matter 
particle (and its siblings) lives in a hidden sector, and inter- 
acts with the visible sector only through heavier particles. In 
this work we have studied several aspects of such a scenario 
in a model independent way (as far as possible). In particu- 
lar, Fermi-LAT is currently probing annihilations into gamma 
ray lines that are 1 or 2 orders of magnitude below the freeze- 
out cross section. If a line is observed in a near future, then 
the hypothesis of annihilation into gluon pairs could naturally 
explain the observed relic dark matter density. This holds in 
particular for the tentative 130 GeV y line claimed in Ref. J6). 
If DM annihilates dominantly into gluon pairs, an inevitable 
by-product is a large flux of anti-protons, which is constrained 
by PAMELA. The upper bounds on the annihilation cross 



section of DM into gluon pairs depend only on the propa- 
gation setup, and not on specific particle model considera- 
tions. Using a standard semi-analytical approach, we give 
in Fig|4]the bounds on (<Tv) gg for three sets of CR propaga- 
tion parameters (MIN, MED and MAX) and two DM pro- 
files (NFW and Einasto, but the dependence on the profile is 
mild). The bounds depend sensitively on the choice of prop- 
agation parameters: a thermal relic candidate is allowed for 
m DM > 180 GeV in the MAX setup, m DM > 80 GeV for 
MED, and m DM > 16 GeV for MIN. More generally, and inde- 
pendently of the relation with thermal freeze-out, but assum- 
ing specific values for the gluon-to-gamma emission ratios, 
Eq.([T]l, the anti-protons data set upper bounds on the annihi- 
lation cross section into y, see Fig. 5. 

Annihilation into gluons also opens the possibility for direct 
detection. The connection between direct detection rates and 
the gluon/photon cross sections is more model dependent, as 
it depends on the topology of the process and the nature of the 
DM candidate. For y-ray lines through a quartic interaction, 
the direct detection rate is in one-to-one correspondence with 
(av) gg . For a scalar DM candidate with WMAP abundance, 
XenonlOO excludes candidates below m^M ,$• 130 GeV, see 
Fig|6] In particular, the predicted rate for the 130 GeV line 
is at the edge of the current exclusion limits. For the other 
topologies, the connection is looser, and they may lead to ob- 



servable direct detection signals, see Section III 



A DM candidate that couples to gluons is constrained by col- 
lider searches, through jets + missing energy signals. Again, 
the correspondence is most straightforward for quartic effec- 
tive operators. Using the results of Ref.[36|, we put in per- 
spective the current and forthcoming constraints from collid- 
ers. As far as we know, the impact of the current data on the 
operators of Eqs.(8]9i has not yet been shown by the LHC 
experiments. It would be interesting to do so. Even if the cur- 
rent bounds are weaker than the current direct detection ex- 
clusion limits for a scalar candidate, they provide independent 
constraints. In the case of the fermionic candidate, collider 
data should set the best current constraints. As for the other 
topologies, we have discussed to which extent they relate to 
the constraints on the quartic interactions. Finally, we have 
put forward the possible production of the intermediate scalar 
(or pseudo-scalar) particle through gluon-gluon fusion, which 
may decay into two photons, a process reminiscent of (light) 
Higgs searches. 

Our results show that it could be interesting to further con- 
strain such a process using available and forthcoming LHC 
data. One last potentially interesting signal (work in progress) 
is the capture of DM in the Sun, and its annihilation into a 
hard spectrum of neutrinos. As emphasized by many authors, 
annihilation into a pair of photons should also annihilate into 
yZ. The latter may decay into a pair of neutrinos, with energy 
E v ~ moni/2 (for a Z-line), oiE v ~ hjom/4 (for a narrow box- 
shaped Z spectrum). 
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